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Abstract. We present N-body cosmological numerical simulations including a primordial non-Gaussianity in the density 
fluctuation field quantified by the non-linear parameter ffj^. We have used MPGRAFIC code to produce initial conditions and 
the Adaptive Mesh Refinement (AMR) code RAMSES to evolve the large scale structure formation. We estimated the higher 
order momenta of the initial distribution of density fluctuations, investigated the redshift evolution of the non-Unear power 
spectrum and estimated the discrepancy introduced by the primordial non-Gaussianity in the non-Hnear power spectrum. 
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INTRODUCTION 

The most stringent constraints on a primordial non-Gaussianity come from the temperature anisotropy measurement 
by WMAP 5-year observations [11], which show a small effect of non-Gaussianity on large scales. On the other hand, 
the probe of a positive non-Gaussianity on smaller scales seems to have been detected recently in WMAP 3 -year 
observations [18]. Extragalactic dust emission can be also associated with primordial deviations from Gaussianity at 
sub-degree scales [15]. 

From the theoretical point of view, the assumption that primordial density fluctuations generating large scale 
structure formation can be modeled by a Gaussian distribution is supported by the inflationary models. However, 
distinct variants of the inflationary models incorporate a primordial non-Gaussianity, as shown by [1]. 

We present a series of N-body numerical simulations of large scale structure formation in order to investigate the 
effect of a primordial non-Gaussianity on the probability distribution function and the two-point correlation function 
of the density fluctuation field of the dark matter distribution. The non-Gaussianity is quantified by the usual parameter 
/nl, which modifies the gravitational potential at a primordial epoch. Higher order modifications of the gravitational 
potential have been recently considered as well [6]. 

NON-GAUSSIAN MODEL 

For a large class of models of generation of the initial seeds for structure formation, including single-field [12] and 
multi-field inflation [2], or the curvaton [9], the non-Gaussianity of the initial density fluctuations can be modeled 
through a quadratic term in the Bardeen's gauge-invariant potential ^, namely [10] 

^(X) = a.L(x) +/nl - i^lix))] , (1) 

where <I>l is a Gaussian random field, with (4>l(x)) = 0, and the specific value of the dimensionless non-linearity 
parameter /nl depends on the assumed scenario. This form of non-Gaussianity can be obtained from a truncated 
expansion of the effective inflaton potential [16]. From the latest WMAP5 temperature data [11] 

- 15 1 < /nl < +253 with 95% CL (2) 

With this convention positive (negative) /nl corresponds to positive (negative) skewness of the probability distribution 
function of density fluctuations. The Gaussianity tests show that the primordial fluctuations are Gaussian to the 0.1% 
level. 

Since the modification introduced in (1) is of second order in the potential, the Unear two-point correlation 
function is not modified for any value of /nl in the chosen interval {fyil < 5,000. In fact, by order of magnitude 
3>(x) < 10~^ — 10~^; so even for /nl ~ 10^ the non-Gaussian correction gives a contribution of the order of 1%. 



It should be noted that distinct reaUzations of primordial non-Gaussianity in the density fluctuation field has been 
proposed, see e.g. [4]. 



N-BODY SIMULATIONS 

The procedure of generation of primordial non-Gaussian initial condition has been implemented in the MPI pa- 
rallelized numerical code MPGRAFIC [14]. The N-body simulations have been performed using the AMR parallelized 
numerical code RAMSES-3.0 [17]. We have followed only the evolution of the coUisionless dark matter particles. We 
performed simulations in computational boxes of physical sizes {I62h~^ Mpc)^, {300h~^ Mpc)^, {324h~^ Mpc)^, 
(500/i^' Mpc)^ with a grid resolution of 256^. The series of simulations presented here aims to detect the non-linear 
scales effect of the primordial non-Gaussianity. 



COSMOLOGICAL PARAIMETERS 

We have adopted the following set of cosmological parameters, within the framework of the so-called "concordance" 
A-CDM model, compatible with WMAP 5-year observations [11]: Q.b = 0.044, Q.^ = 0.26, Q.a = 0.74, H = 72 km 
s^^Mpc^^ The spectral index of the initial power spectrum is n = 0.951. The rms of density fluctuations in spheres 
of radius Rq, a^{Ro), is defined as 

a^{Ro) = J d^kP{k;z)w\kRo) (3) 

in absence of non-linear perturbations, i.e. P{k; z) is the linear power spectrum. Here {kRo) is an appropriate window 
function. At the scale Rq = Sh'^ Mpc, we imposed c7^(^o) ~ 0.792 [11]. 



INITIAL CONDITIONS 

We first briefly review the procedure of generation of 3D Gaussian initial conditions [3, 14]. In the Gaussian case, 
the 3D initial density and displacement fields are completely defined by the two-point correlation function, or power 
spectrum, P{k) of the density fluctuation field 5(x), defined as /'(A:)5D(k — k') = (5(k)*5(k')), where 5d(x) is the 
Dirac delta function and 

5ik) = j^Jdxe-''^5{x) (4) 

The 3D density fluctuation field is generated as convolution of a Gaussian white noise n(x) with the CDM transfer 
function T{k) = ^JP{k). In the position space, the Gaussian field n(x) is generated with the standard technique of the 
Gaussian white noise, using a uniform phase and an amplitude extracted from the Rayleigh distribution. The field n(x) 
has zero average value and unitary variance, i.e. (|n(x)p) = 1, or, equivalently, constant power spectrum: (|n(k)p) = 1. 
As a second step, the Gaussian white noise n(x) is convolved with a CDM transfer function T{x) to obtain the density 
fluctuation and displacement fields according to the Zel'dovich approximation [19], i.e. 

5(x) = j dke*''«(k)r(yfc), Sj{x) = J dke*''"^'^^P^%. (5) 

with J = 1,2,3. Therefore, the primordial fields 5(x) and Sj{x) are made up by two decoupled contributions: the field 
n{x) represents the statistics at early epoch, which is commonly assumed to be Gaussian, providing the phase of the 
fluctuation. On the other hand, the unknown physics of plasma processes leading to the CMB epoch is contained in the 
transfer function T{k), completely determined by the two-point correlation function, for which the parametrization of 
[7, 14] is used. The T{k) provides the amplitude of the fluctuation and depends on the cosmological model. 

In order to disentangle primordial non-Gaussianity from non-Gaussianity generated by the non-linear evolution of 
structures under the gravitational potential, the statistics has been modified at early times, i.e. before the convolution 
with the transfer function, by using the (1). Once the Gaussian white noise «(x) is generated, the primordial den- 
sity fluctuations field in the Fourier space 5^(k) = n(k) is related through the Poisson equation to the primordial 
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FIGURE 1. PDF for different values of fyi at redshift Zstart ~ 77. 



gravitational potential <I>''(x): 

^2^P(k) = ^5'^(k) (6) 

where £l„^^, i.e. the matter density in units of ciitical density, and = c/Hq are evaluated at the present epoch. The 
potential <I>''(k) is first transformed back to the position space, and converted in a non-Gaussian potential according to 
the prescription in (1). The new gravitational potential <I>^p(x) coiTesponding to a non-Gaussian primordial statistics 
is converted back to the Fourier space in which the new density fluctuation and displacement fields are computed. 
The initial condition for the N-body simulation is obtained by using the new non-Gaussian gravitational potential to 
perturb the initially uniform density distribution according to the Zel'dovich approximation. 
The initial redshift, as input for the code RAMSES, is Zstan ~ 77. 

In Fig. 1 the probability distribution function PDF(5) is shown for different values of fyi. The PDF(5) of the 
fluctuation density field 5 at a given scale R is computed by sampling the simulation box at a given redshift with a 
number of independent spheres of radius R, fully covering the simulation box. In Fig. 1 the PDF with comoving sphere 
radius R ~ 8/?^' Mpc is shown in a simulation box of size L = 162 h^^ Mpc and initial grid resolution 256^, with a 
cell grid comoving size of AL = 0.63/;^' Mpc. 

Significant discrepancies with respect to the Gaussian case are found in the tail of the PDF only for values of fyi 
beyond the constraints inferred from WMAP 5-years [11], at variance with [8, 13]. 

The second and third order momenta of the 5 distribution have been estimated at the initial redshift Zstart ^ 77, 
corresponding to the initial condition existing prior the gravitational evolution under the hydrodynamics equations. 
We find that at Zstart. for fyi = 5000 and /atl = —5000 the third order momentum is 8^ = 4.61 x 10"** and 8^ = 
—3.93 X 10"^, respectively. 

As a check the numerical reliability of our computation, the third momentum of the non-Gaussian density fluctuation 
8 can be easily estimated, as pointed out by [5]. From the Laplacian of (1), V^O = V^Ol + 2/nl [| VOlP + 4>lV^<I>l] , 
by using the Poisson equation, the last relation implies (5^) — 6fNL{8l'i>L), at the lowest order in the product /nl'^'l- 
This relation has been verified with good accuracy. 
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FIGURE 2. Left panel: Redshift evolution of the dimensionless power spectrum l^{k) toward z = 0, for different values of ff^i. 
Right panel: Comparison at z = of the power spectrum of matter density P{k) in the Gaussian and non-Gaussian case. We find a 
discrepancy of the order of 15% (2 - 3%) at = 0.7(/i Mpc" ' ), for fNL = ±5000 (±750). 



DEFINITION OF /nl 

As already pointed out by several authors, see e.g. [11, 14], the definition of /ml depends on the cosmic epoch (i.e. on 
the value of the FRW scale factor a) at which (1) is applied, because both the potentials Gaussian and non-Gaussian, 
<I>L and <I> respectively, evolve in time with the factor g{a) = D{a)/a, where D{a) is the linear growth factor of density 
fluctuations. Therefore the relation between the fyi^ introduced at z — > oo and fyi^ introduced at present epoch is 
/nl^ = fNLQg{0)/g{°°), where g{°°)/g{0) ^ 1.34. In the simulations presented here, the gravitational potential is 
modified directly at primordial epoch z — > 0°, before the convolution with the CDM transfer function T{k). This choice 
allows for a direct modification of the primordial statistics (see Sect, above). Therefore, no linear extrapolation back 
to Zstart of the gravitational potential is peiformed. The parameter /a?/, introduced in (1) affects only the statistics of the 
matter distribution and does not affect the physics pre-CMB epoch which is included by the further convolution with 
the transfer function. 



RESULTS 

The formalism summarized above allowed us to numerically explore the imprint of the primordial non-Gaussianity 
of fNL-type on the dimensionless power spectrum of the CDM density field A'^{k) ~ P[k)k^ / {2n^) at z = 0, see left 
panel of Fig. 2. The computational box has a size of 162 h^^ Mpc, with an initial coarse graining of 256^ cells. Higher 
resolution simulations are currently underway. The chosen resolution allowed us to resolve the small scale contribution 
of the non-Gaussian perturbation to the structure formation. 

At z > 5 the discrepancy between the Gaussian and non-Gaussian case decreases because non-linear effects due to 
primordial non-Gaussianity on the scale of interest become negligible at large redshifts. 

In the right panel of Fig. 2 the ratio of power spectra P{k; fyi) / P{k; fyi — 0) is extracted from the simulations 
at redshift z = 0. The matter power spectrum of non-Gaussian models appears to deviate by a few per cent at 
k = O.lhMpc^^ , with /^f^ = ±5000. Moreover, deviations from Gaussian case are enhanced at non-linear scale. 

CONCLUSIONS AND PERSPECTIVES 

The formalism of /A^i-type allowed us to numerically explore the imprint of the primordial non-Gaussianity on 
a number of observables quantities, i.e. the initial probability distribution function of the density fluctuation field 
PDF {5) and the non-linear power spectrum at present epoch. 



The amplitude of departure of power spectrum from the Gaussian case could be hardly detectable, also due to a 
progressive coalescence of power spectra as z decreases. 

The detectabihty of third order momentum of 5, regardeless the very small amphtude, can lead to observable 
differences at the present epoch in the case of initial predominance of sub-density regions. 

A series of simulations with higher space resolution will be performed to analyse the mass function and the bias 
of dark matter haloes within the same model for the primordial non-Gaussianity at the leading perturbative order in 
the gravitational potential, with parameter fyi. On the other end some simulation tests involving the successive term 
in the Taylor expansion of the potential have been carried out, with the parameter g^i [11], leading to not-observable 
deviations from Gaussian case in the PDF and P{k;z = 0) up to ^^vl =10^. 
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